Protein Design Using Continuous Rotamers

Optimizing amino acid conformation and identity is a central problem in computational protein design. Protein design algorithms must allow realistic protein flexibility to occur during this optimization, or they may fail to find the best sequence with the lowest energy. Most design algorithms implement side-chain flexibility by allowing the side chains to move between a small set of discrete, low-energy states, which we call rigid rotamers. In this work we show that allowing continuous side-chain flexibility (which we call continuous rotamers) greatly improves protein flexibility modeling. We present a large-scale study that compares the sequences and best energy conformations in 69 protein-core redesigns using a rigid-rotamer model versus a continuous-rotamer model. We show that in nearly all of our redesigns the sequence found by the continuous-rotamer model is different and has a lower energy than the one found by the rigid-rotamer model. Moreover, the sequences found by the continuous-rotamer model are more similar to the native sequences. We then show that the seemingly easy solution of sampling more rigid rotamers within the continuous region is not a practical alternative to a continuous-rotamer model: at computationally feasible resolutions, using more rigid rotamers was never better than a continuous-rotamer model and almost always resulted in higher energies. Finally, we present a new protein design algorithm based on the dead-end elimination (DEE) algorithm, which we call iMinDEE, that makes the use of continuous rotamers feasible in larger systems. iMinDEE guarantees finding the optimal answer while pruning the search space with close to the same efficiency of DEE. Availability: Software is available under the Lesser GNU Public License v3. Contact the authors for source code

Structural Elements Regulating Amyloidogenesis: A Cholinesterase Model System

Polymerization into amyloid fibrils is a crucial step in the pathogenesis of neurodegenerative syndromes. Amyloid assembly is governed by properties of the sequence backbone and specific side-chain interactions, since fibrils from unrelated sequences possess similar structures and morphologies. Therefore, characterization of the structural determinants driving amyloid aggregation is of fundamental importance. We investigated the forces involved in the amyloid assembly of a model peptide derived from the oligomerization domain of acetylcholinesterase (AChE), AChE586-599, through the effect of single point mutations on β-sheet propensity, conformation, fibrilization, surfactant activity, oligomerization and fibril morphology. AChE586-599 was chosen due to its fibrilization tractability and AChE involvement in Alzheimer's disease. The results revealed how specific regions and residues can control AChE586-599 assembly. Hydrophobic and/or aromatic residues were crucial for maintaining a high β-strand propensity, for the conformational transition to β-sheet, and for the first stage of aggregation. We also demonstrated that positively charged side-chains might be involved in electrostatic interactions, which could control the transition to β-sheet, the oligomerization and assembly stability. Further interactions were also found to participate in the assembly. We showed that some residues were important for AChE586-599 surfactant activity and that amyloid assembly might preferentially occur at an air-water interface. Consistently with the experimental observations and assembly models for other amyloid systems, we propose a model for AChE586-599 assembly in which a steric-zipper formed through specific interactions (hydrophobic, electrostatic, cation-π, SH-aromatic, metal chelation and polar-polar) would maintain the β-sheets together. We also propose that the stacking between the strands in the β-sheets along the fiber axis could be stabilized through π-π interactions and metal chelation. The dissection of the specific molecular recognition driving AChE586-599 amyloid assembly has provided further knowledge on such poorly understood and complicated process, which could be applied to protein folding and the targeting of amyloid diseases

Pan-cancer analysis of whole genomes

Cancer is driven by genetic change, and the advent of massively parallel sequencing has enabled systematic documentation of this variation at the whole-genome scale(1-3). Here we report the integrative analysis of 2,658 whole-cancer genomes and their matching normal tissues across 38 tumour types from the Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium of the International Cancer Genome Consortium (ICGC) and The Cancer Genome Atlas (TCGA). We describe the generation of the PCAWG resource, facilitated by international data sharing using compute clouds. On average, cancer genomes contained 4-5 driver mutations when combining coding and non-coding genomic elements; however, in around 5% of cases no drivers were identified, suggesting that cancer driver discovery is not yet complete. Chromothripsis, in which many clustered structural variants arise in a single catastrophic event, is frequently an early event in tumour evolution; in acral melanoma, for example, these events precede most somatic point mutations and affect several cancer-associated genes simultaneously. Cancers with abnormal telomere maintenance often originate from tissues with low replicative activity and show several mechanisms of preventing telomere attrition to critical levels. Common and rare germline variants affect patterns of somatic mutation, including point mutations, structural variants and somatic retrotransposition. A collection of papers from the PCAWG Consortium describes non-coding mutations that drive cancer beyond those in the TERT promoter(4); identifies new signatures of mutational processes that cause base substitutions, small insertions and deletions and structural variation(5,6); analyses timings and patterns of tumour evolution(7); describes the diverse transcriptional consequences of somatic mutation on splicing, expression levels, fusion genes and promoter activity(8,9); and evaluates a range of more-specialized features of cancer genomes(8,10-18).Peer reviewe

Anatomy of the Hamstrings

This chapter will provide the anatomical foundation for the content to come in later portions of this book. It will begin with an overview of the proximal insertion sites of the muscles that comprise the hamstring group. The proximal tendons and musculotendinous junctions (MTJs) of semimembranosus, semitendinosus and the long and short heads of biceps femoris long head will then be described, highlighting the differences in structure between each of the muscles. The distinct architectural characteristics of each muscle belly (e.g. size, fascicle orientation within and between muscles) will be outlined, followed by the structure of the distal tendons and MTJs. Finally, a summary is provided of the neurovascular supply of the hamstrings